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ABSTRACT 


The scope of this thesis is to investigate the vertical and horizontal plane motions 
of surface ships in close proximity towing in irregular waves. Strip theory calculations 
have been utilized in order to predict the hydrodynamic coefficients and wave exciting 
forces and moments in sway and yaw, heave and pitch. The appropriate matching 
conditions between the two ships are provided in terms of the resistance-speed 
characteristics of the leading ship. The two-parameter Bretschneider spectrum with a 
cosine-squared spreading function is used to model the sea state environment. An 
extensive set of parametric studies is presented in a wide variety of developing and 
decaying sea states. 
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I. INTRODUCTION 


The ocean waveform is an undulating process that transmits energy thousands of 
miles both efficiently and effectively. A twelve-second periodic wave group will take 
forty-eight hours to cross 1,000 miles of deep, open ocean with little loss of energy. 
When its power is unleashed on a rigid body, it becomes one of the greatest fears of all 
mariners and maritime insurance companies alike. It is the effort of the naval architect to 
design ships that can surpass a broad spectrum of waveforms and the energy associated 
with them. 

A ship advancing at a relatively constant forward speed with an arbitrary heading 
in a train of regular waves will move in six degrees-of-freedom (three translational: 
surge, heave, sway; three rotational: pitch, roll, yaw). In order to be able to compute all 
responses of the vessels to random seaways, it is necessary to deal with the complete 
motions of the vessel in all six degrees-of-freedom, with special attention to the important 
couplings of heave and pitch, sway and yaw. For such a rigid body, six nonlinear 
equations of motion, with six unknowns, must be set up and solved simultaneously. In 
the case of a vessel of port/starboard symmetry and restrained about this axis of 
symmetry, the six non-linear equations can be reduced to two sets of three linear 
equations. 

Specifically, it is the effort of this thesis to present a linear theory of ship motions. 
While most vessel responses are nonlinear to some extent, however, where nonlinearities 
are small a linear theory will yield reasonable predictions. Experimental and theoretical 
investigations have repeatedly shown that a linear theory analysis gives excellent 
predictions over a wide variety of sea states and vessel platforms (Beck). 
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A. MOTIVATION 

Since the advent of steam, the level of engineering necessary to keep ocean-going 
vessels afloat and righting has had to become more and more complex. The economic 
demand of trading nations, and the warship designs of international navies continually 
pushes the complexity level to yet another extreme. Since the rapid developments in 
hydrodynamic theory of the 1950s, the advances in oceanography and computer 
technology have enveloped the marine engineers task in pursuing advanced designs in 
high-speed hull forms, and unusual platforms such as SWATH vessels. 

Close proximity towing of surface ships is a matter of interest to the U.S. Navy 
and the Office of Naval Research. Several possible applications have been suggested 
including the SLICE/KAIMALINO connection (Nash) and the SEA LANCE 
configuration (TSSE). Of primary concern in this dissertation is the evaluation of motions 
in six degrees-of-freedom of the two ships when in close proximity towing operations. 
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II. WAVE SPECTRA 


A. EARLY WAVE THEORY 

Franz Gerstner of Czechoslovakia officially documented the first, rather primitive 
wave theory in 1802. He described how water particles move in oscillatory motions 
within the wave. He further described how the water particles in the crest move in the 
direction of wave advance and those in the trough move in the opposite direction. Such 
observations had been made centuries earlier; it was his discovery, however, that 
described the circular path of the water particles to have a diameter equal to the height of 
the passing wave (path diameter proportional to water depth). He essentially described a 
trochoidal wave (Figure 1). This discovery bridged the gap between the casual seafarer’s 
observation to that of mathematical derivation and hence the theoretical beginning of 
hydrodynamics. 

Following the theoretical 
came the experimental. The first 
experimentalists developed 

controllable wave tanks with the 
capability to reproduce waveforms 
over and over, and whereby careful 
measurements could be taken. 

Then in the 1950’s a new breed of 
experiment was developed 
harnessing modern computations 
and statistical applications; and, 
hence, the numerical wave tank was 
bom. 
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GEOMETRY OF TROCHOID 


Figure 1. Trochoidal Wave (After Beck) 
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B. WAVE SPECTRA ANALYSIS 


The understanding of vessel response at sea, and the ability to predict their 
behavior begins with the study of the nature of ocean waves in which the vessels will 
operate. Consequently, in order to understand the behavior of ocean waves, one must 
first understand their origin. Waves are a natural phenomenon that occurs between the 
interface of any two fluids of different density. While any kind of disturbance in the 
ocean will generate waves of sparing severity, there exist four primary sources: wind, 
earthquakes, landslides, and the gravitational attraction of the moon and sun. Of the four, 
the greatest of importance is the wind-generated waves. The size and variety of these 
wind waves depends on many factors, namely, the velocity of the wind, the distance it 
blows across the water, the time duration that it blows, and the depth of water in which 
the wind transverses (Figure 2). The character and signature of the wave is quickly lost 
as the wind energy is transferred to tidal energy in this highly stochastic ocean media. 
Once the winds have transferred their energy to the ocean surface, the waves begin their 
natural procession from the generation area. The original wind waves eventually decay, 
evidenced by the crests as they become lower, and more rounded and symmetrical. 
These crests then begin to move in groups of similar period and wave height, and their 
form becomes more and more sinusoidal. This origination theory becomes now the basis 
for our analytic discussion: sinusoidal wave forms moving as steady-state, constant 
parameter wave “trains” under the control of gravity and inertia; and their subsequent 
influence on sea-going vessels in-tow. 
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Wave Frequency, rad/sec 

Figure 3. Sea State Four (Varying Period) (From: McCreight) 

Waves are classified according to their period, which can range from a fraction of 
a second to more than ten thousand seconds. They can range in size from a ripple to 
great storm waves and tides whose wavelength can reach half the circumference of the 
earth (Bascom). Essentially we are studying a medium whereby the wave spectrum 
ranges from waves so small they can barely be seen to waves so long that they are not 
even noticed. For these extreme conditions our vessel design is safe. It is a narrow band 
of wave spectrum that can cause greatest damage to seagoing vessels and is of utmost 
concern to the naval architect and our research. 

The energy in the ocean is distributed among several distinct groups of waves 
each with a defined range of periods (Figure 3). All waves can be classified as gravity 
waves since once they are started; gravity is the driving force that keeps these sea crests 
moving in a never-ending attempt to restore the original flat sea surface. The potential 
energy of each wave is due in part to the waves profile which consists of two sides of the 
equilibrium position: the crest that rises above datum sea line, and the trough that 
extends below. A wave group moves perpetually attempting to overtake the trough ahead 
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and restore equilibrium, suffering of course from the wave ahead with the same intent. 
For swell, the wave’s potential and kinetic energies are equal. The kinetic energy 
possessed by a wave is due to the imparted motion of the water particles, while the 
remainder of the waves energy is the potential energy due to the elevation of the center of 
gravity of the mass of water in the crest above mean sea level (Figure 4). 



Another classification of waves, although when fully developed act like gravity 
waves, are called capillary waves (short crested seas). They are a special waveform 
whose driving force is not gravity, but rather the surface tension at the water/wind 
interface. The capillary force is stronger than gravity, and given the open ocean 
condition, thrives in producing the most abundant kind of wind-generated wave. 
Capillary waves give rise to the development of ripples (incident wind drag), which are at 
the lower end of the wave spectrum, but none-the-less lead to the growth and overall 
combined effect of fully developed seas (wave energy absorption proportional to wind 
velocity). 

Open ocean waves do not have the normal and precise properties as those 
generated in a wave generator. The height of the crests and the depths of the troughs are 
highly irregular. The theoretical terminology of a discernable wave period and wave 
velocity (celerity) is lost when taken to the real and observable, highly stochastic and 
irregular environmental oceanic swirl. For these very reasons is why the physics of 
ocean waves must rely on statistical methods in order to describe the properties of waves 
(Appendix B). 

Wave heights, while apparently just as irregular, can be referenced to a datum line 
also referred to as mean sea level. Using mean sea level as a baseline, a sea state is 
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simply the result of superimposing a number of sinusoidal wave trains one on top of the 
other (Figure 5). 

Each layer represents a series of regular sine waves, each having their own 
distinct wave height, wavelength, and direction; and coincides with the theoretical 
notions exactly. The summation of the wave crests and troughs of varying amplitude 
tend to cancel each other out, and conversely, amplify the local maximum or minimum 
depending on coincidence of discrete wave points throughout the ocean surface. Again, 
as these waves transverse each other, the cancellation and amplification will subside 
quickly as these respective wave groups transverse their own direction. Accordingly, the 
greater the number of wave layers the more random the sea surface. 

The seemingly random nature of ocean waves can be better characterized by its 
energy spectra. This methodology assigns a value to the square of the wave height for 
each frequency and direction. Then, using statistical methods, the wave spectrum can be 
analyzed to reveal how the wave energy is distributed among various wave periods. 



Figure 5. Wave Fourier Composition (From: Ochi) 
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The essential feature of these mathematical techniques/models is the concept of a 
spectrum defining the distribution of energy among the different hypothetical regular 
wave components having various frequencies (wave lengths) and directions. It has been 
formulated that the irregular motions of a ship in a seaway can be better described as the 
linear superposition of the responses of the same ship to all the wave components of the 
seaway. For our purposes, the linear theory assumption, as applied to the SWATH 
vessels, implies the vessels’ response is described by the same statistical properties of a 
random (stochastic) surface wave environment. Given these assumptions, these surface 
waves can be characterized as having a Gaussian, or normal, probability distribution 
about some mean (zero amplitude, calm seas), under short-term statistically stationary 
conditions (Figure 4). This statistical technique will greatly simplify the numerical 
involvement of our study with regards to the application of statistics, probability theory 
and Fourier analysis techniques. 

During the energy transfer between wind and water, there are several important 
wave interactions and wave-breaking processes that affect the dispersion and propagation 
of waves from the storm area. Particularly for the case of short crested, small amplitude 
waves, the principle of linear superposition applies. For example, if Xi(xi,yi,ti) and 
^ 2 (x 2 ,y 2 ,t 2 ) are two wave systems, then Xi(xi,yi,ti) + X 2 (x 2 ,y 2 ,t 2 ) is also a wave system. 
Inherent to this assumption, the two systems can move through each other without 
damping the other whereby all wave energy is conserved linearly. Of equal conceptual 
concern is wave celerity, or velocity, which is wholly a function of wavelength (shorter 
waves travel slower than longer waves). Referring to Figure 5, it can be shown that any 
wave system can be divided into a sum of component regular waves of various 
wavelength, amplitude and direction using Fourier Integral techniques. Since water is 
considered to be incompressible, the average value of vertical displacement at any instant 
for a sinusoidal, or regular, wave period is zero. While the mean will give us little 
statistical information, the variance, however, is a positive value that will directly relay 
the severity of the sea. A fundamental theory of statistics states that the variance of the 
entire system is determined by measuring the amount of variability in the distribution of 
the independent wave groups themselves. In general, the variance of a continuous 
function can be expressed by the following: 
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(SOO 2 ) = 2 {t)dt (i) 

And in the case of simple harmonic sine waves, the variance of wave elevation of 
a single cycle is equal to one-half the square of the wave amplitude. 

The slight errors involved in the assumed sinusoidal wave shape of linear 
harmonic wave theory are standard practice, and can be considered reasonable 
assumptions with little consequence. The errors only become significant when dealing 
particularly in the wave breaking phenomena, and when the geometry of the wave crests 
becomes unsteady. The correct mathematical analysis of linear short-crested irregular 
waves is formable with little error using short-term stochastic models of ocean waves. 

Wave Properties (Two-Dimensional Waveform): 

C (x,y, t) = C, cos[^(vcos P + ysin p )-a>7+e j 


£=phase angle 

P=wave direction (global coordinates) 
to=wave frequency 
Wave Number: k - 2n / L w = oo 2 / g 


V 2 

Wave Length: L w = 271 — 1 — 

§ 


Wave Period: T = 


r 2nL ^ 


8 J 


Surface Profile: Co = C cos k(x-V c t) 


Velocity Potential: (f> = ~C V c e :z sin k{x-V c t) 


Wave Variance: 


Wave Celerity: 



V 


K 

T 


vf 

u 


(Appendix B) 


9 



Continuing with the short-term model, the total wave system can be shown as a 
summation of many, if not infinite, independent wave components: 

C(0=LC~cos(-GV + e ; ) (2) 

i 


This system can also be represented by a point spectrum, also known as a 
variance spectrum, S(co) (the wave is observed only at a point, and none with respect to 
wave direction). At any particular wave frequency, CO;, the variance of all the wave 
components within a small finite frequency band, 8co, centered at 0); is denoted by: 

(C 2 ) = S'(co,. )8 co (3) 

following: 

(v)=y 2 v (4) 

substituting: 

y 2 t; 2 = S{ co,.)8co (5) 

Integrating the RHS of the equation will resolve a good approximation of the 
wave system’s total energy (Figure 3). 


E =Js (oo)Jco (6) 

0 


In the derivation of energy spectra in short crested seas, Figure 6 depicts only 


fourteen wave groups, a multi-directional short crested sea requires many more (73 was 



SCSi E OF TIME OH Otsr Mice 

Figure 6. Variance Spectrum of Waves (From: Beck) 
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Figure 6, depicts several distinct sine curves of random phase. Numerical 
summing these will result in the irregular wave pattern in which we are interested. Figure 
5 demonstrates this idea quite clearly. 

While this seems a satisfactory depiction of the development of ocean wave 
spectra, and mathematically the variance, E, obtained from a point spectrum is a decent 
measure of the sea severity, it can be further shown that this measurement is still 
insufficient. A more complete characterization of the seaway must include wave 
component direction, also known as a directional spectrum. From the two-dimensional 
case to the three-dimensional case our generalized wave equation becomes: 

C (x,y,t) = IE? ij ;Cos[k,(xcos p ; + ysinP^-cqr+e^] (7) 

i j 

and 

(c (o 2 ) ^ e= y 2 ££57=j' £” s«o ,p xprfm (8) 

« j 

Wave energy for a point has an angular distribution as well as a distribution over 
a range of frequencies. This angular distribution of wave energy is termed directional 
spreading and overall greatly increases the accuracy of our results when designing for 
short crested sea behavior. The directional spectrum thus defines the distribution of 
energy with even greater resolve defined by the frequency and direction of the variances 
of the individual wave system components (each having a unique combination of 
frequency and direction, and random phase angle). Furthermore, for our incorporated 
spreading function, the directional spectrum can be subdivided into the following: 

S(g>,P) = S«o)M(P) (9) 

where: 

M(P) = —cos 2 P (10) 

71 

While the frequency and direction are accounted for, what is critical to our study 
is defining an idealized random seaway. We have already discussed the constructs of 
wave spectra from recorded data, but how can we reconstruct such data given certain 
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environmental parameters? This technique is particular to the design process, where we 
are only interested in recreating wave energy spectroscopy. This approach allows the 
designer to modify such environmental parameters incrementally to discern ship response 
to specific wave parameters. 

The Bretschneider Spectrum was selected and can be presented as follows: 

l S(co) = ^^//V 125(<n '" / “ )4 (11) 

4 co 


(0 m 


2 TZ 

T 

m 


( 12 ) 


where: 


T m =input modal wave period 

H s =input significant wave height defined as the average 
one-third highest wave heights (H 1 /. 3 ) 

This is a two-parameter wave spectrum and is a widely accepted design parameter 
for marine systems. Using this infamous two-parameter spectral formulation it is 
possible to generate a myriad of wave groups of varying severity and size composition. 
This spectral density function can quickly illustrate the magnitude of wave energy. The 
most commonly used definition of sea severity is significant wave height, which can be 
deduced from the wave spectrum as being equal to four times the square-root of the area 
under the curve (Figure 3). 

Continuing with the linear theory assumption, the spectral density of any response 
can be found by multiplying the incident wave spectrum by the square of the response 
amplitude operator (RAO) of the desired response. In control theory, for example, the 
RAO is often called the transfer function of the linear system function. For any given 
frequency, the RAO is the amplitude and phase of the desired response to regular incident 
waves. Also known as the motion transfer function, the response amplitude operator 
(RAO) maps the complex response of a vessel to a seaway or input spectrum as a 


12 



function of frequency. The spectrum of the particular response that the RAO has been 
computed for is calculated from: 

S R (a) - |RAO(a))| 2 S(oo,p ) (23) 

In this equation, S R ( 03 ) is the response of the vessel to the input sea spectrum for 

a given frequency. Ship response then can be calculated for a given wave frequency and 
significant wave height. For example, the complex absolute motions predicted in regular 
wave modeling (? s , ?k) are converted into RAO’s for absolution motion as 

RAO(^ sk ) = abs(^ sk ) (24) 

The response spectrum for absolute motion then becomes: 

S R -^ s , k m = \ abs ^s, k ) | 2 S(co,{3) (25) 

Finally, all motion analysis speed polar plots are resolved in RMS values: 

(RMS) 2 = J" Js r (26) 

P CO 
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ID. SHIP RESPONSE 


The longitudinal motions of a symmetrical ship with symmetrical moorings in 
regular waves can be considered distinctly from the horizontal degrees of motion. The 
principle modes are the heave and pitch degrees-of-freedom (modes 3 and 5, 
respectively). In particular, the cross coupling between these two modes is critical to the 
success of this model. To further this discussion, it has also been shown that the surge 
degree of freedom (mode 1), given the SWATH vessels’ architecture, has only a minor 
contribution and can be neglected for the present discussion (refer: Conclusions pp.23) 
For the following derivation, the SLICE/KAIMALINO coupling is rigid, has constant 
forward velocity, and are resolved into the simplified head seas case (waves from the 
bow; (1=180°). It is also assumed that the wave excitation forces are linear and harmonic, 
acting according to the wave encounter frequency (Poppler Shift). Using the linear 
theory assumption, the responses of the vessel will be directly proportional to wave 
amplitude and the frequency of occurrence corresponds to the frequency of encounter. 
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From Newton’s second law of motion, a translational displacement is defined as 
the mass times the acceleration, therefore, for heave we have (center of gravity located at 
the waterline): 

AtT 3 =F 3 (27) 

and for pitch, the rotational modes equal the mass moment of inertia times the angular 
acceleration: 

I 55 % = F 5 (28) 

where A equals displacement, I55 is the mass moment of inertia about the y-axis and F 3 
and F 5 are the total force and moment acting on the body. 

For our discussion (linear theory), the total force and moment consist mainly of 
hydrostatic and hydrodynamic fluid forces: 

F 3 (t)=Fex3(t)+F H 3(t) (29) 
and 

F 5 (t) = F gx 5 (t)+F H 5 (t) (30) 

where F ex is the exciting force due to the waves acting on the restrained ship, and Fh is 
the radiation force due to ship motions in an assumed calm sea, ideal case. Again due to 
a linear theory assumption, the fluid forces can be divided as such. 

— |^£X3 |cOS(CO c t +£3) (31) 

^EX5^ ~ |^Et5| C0S (®f/ £ 5 ) (32) 

where IFex 3 I refers to the amplitude of the heave force and IFrxsI refers to the amplitude 
of the pitching moment. Also called the Froude-Krylov, or Korvin-Kroukovsky, 
excitations (NASH), these excitations represent the integration of the pressure field, over 
the ship wetted surface area, which would have existed in the incident wave system if the 
ship were not present. The diffraction excitation forces and moments, which are also a 
component of IFex 3 I and IFrxsI- are caused by the diffraction or modification of the 


16 



oncoming waveforms. The respective epsilon values are the phase angles between the 
excitation force and sea waves. 

F Hi = [A 33 ((oXi ' 3 + B 33 ( oo)ri 3 + C 33 (co)r | 3 + A 35 (co)ti ' 5 + B 35 (&)r\ 5 + C 35 (co)ri 5 (33) 

F h5 = [A 53 (co)ri ' 3 + 5 53 (to)rj 3 + C 53 (co)r | 3 + A 55 (to)ti ' 5 + 5 55 (co)rj 5 + C 55 (co)ri 5 (34) 

where A,i-terms correspond to added mass in phase with vertical accelerations, the Bjk- 
terms to hydrodynamic damping in phase with vertical velocity, and the Cjk- terms to 
restoring forces (buoyancy effects) in phase with ship displacements. These coefficients 
were developed using a strip theory (a fast and reliable numerical method, capable of 
accommodating a wide range of hull forms). (A 35 (co)ri ' 5 represents the force in the heave 
mode due to an acceleration in the pitch mode.) 

As may be evident to the reader, the coupled equations are similar in form to that 
of a two-degree-of-freedom mass-spring-damper system (added mass+mass, damping, 
and restoring terms on the LHS; excitation forces on the RHS). The dynamics of the two 
systems, however, are far from similar. For the ship system all coefficients and 
excitation forces are all functions of frequency, whereas for the simple mass-spring- 
damper system these values are time independent. While solutions to the vessel coupling 
problem are possible in the time domain, such methods require heavy mathematics 
involving complicated convolution integrals. In order to overcome this difficulty the 
problem need only be solved in the frequency domain. The exciting force, therefore, can 
now be defined in the frequency domain as: 

F ex (t)=F ex e^ (35) 

Defining the motion and its derivatives: 

ii(f) =i]e mJ 

t)(0 = (36) 

f\(t) = -w V”'' 

where rp is the complex response amplitude containing both magnitude and phase of the 
response. 
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Finally, linearizing the equations of motions (Abkowitz, 1969), (note Figure 1): 


Mr\ l + z c r\ 5 ) = F l 
A0i' 2 - z c r| 4 + x c t\ 6 ) = F 2 


AOij-^r \ 5 ) = F 3 

^44^4 — ^46^6 _ A Z( Ti 2 — F 4 

■^5 5*15 A t c fil _ W — ^5 


(37) 


^66^6 7 6 4n 4 + A.r c Ti 2 — F () 


Rewriting the equation of motion in generalized form: 


£[-«);(A„ + A ) + im,B Jk + C^M, = FJ + Ff 


(38) 
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IV. RESULTS 


The following results (Table 1) are presented in speed polar plots (Appendix A). 
As discussed earlier, waves on the bow, or head seas, correspond to 180° on the plots. 
The input parameters as necessitated by the Bretschneider Equation are modal period and 
significant wave height, with the additional input of connection ratio on behalf of vessel 
coupling t> determine connection force and, ultimately, modal excitation. The radial 
coordinates range from zero to twenty knots, and all results show contours of constant 
connection force (RMS values). 

Table 1. Numerical Test Matrix 



1 

Modal Period 

Connection Ratio 


5 

5 

0.1 

0.5 

1.0 


5 

10 

0.1 

0.5 

1.0 

5= 

O) 

5 

15 

0.1 

0.5 

1.0 

'o 

x 

10 

8 

0.1 

0.5 

1.0 

CD 

IQ 

12 

0.1 

0.5 

1.0 

03 

> 

10 

16 

0.1 

0.5 

1.0 

5 

20 

8 

0.1 

0.5 

1.0 

c 

03 

o 

"F 

20 

12 

0.1 

0.5 

1.0 

20 

16 

0.1 

0.5 

1.0 

CD 

C 0 

30 

12 

0.1 

0.5 

1.0 

3Q 

16 

0.1 

0.5 

1.0 


30 

20 

0.1 

0.5 

1.0 


Somewhat illusive to the discussion has been from where are such motions 
observed? Up to this point it has been assumed that the wave spectra has been with 
respect to some fixed referenced point. The absolute displacements, velocities, and 
accelerations have all been assumed to be with respect to an absolute or non- 
rotating/transversing observer. However, if the velocity of the wave groups along with 
the velocity of the vessel is taken into account the results will change significantly. The 
time record, or frequency of encounter, of a moving point relative to a moving fluid 
media and some global fixed inertia reference frame will undoubtedly be significantly 
affected by the Doppler Shift in the component frequencies of the wave system itself. 
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Figure 8. Frequency of Encounter (From: Ochi) 


The location of the ship can be given by: 


x o =X+U a t 

y a = y 


(40) 


Where the component U 0 t is the difference between the body fixed frame and the 
global coordinates. Substituting into the two-dimensional wave equation, the expression 
for the wave field as observed from the moving ship becomes: 

£ (x, y,t ) = £ cos[kxcos p + ky sin p - (to - kU a cos p)t +£ ] (41) 


From Figure 8 the frequency of encounter can then be described as: 


U W ~U 0 cosp 


(42) 


w, = —(U w - U a cos p) (43) 

CD 2 

= ® 0 -—t/ G cosp (44) 

§ 

co e = co(l - a) (44) 
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For the simplified solution: 

a<0: bow waves, co e >co 

0<a0.5: stem waves co e <co. (waves pass the ship quickly) 
maximum co e when a=0.5 

0.5<a<l: stem waves co e <co (waves pass the ship slowly) 
a=l: wave velocity equal to that of the ship. 

a>l: waves come from the stern, but the ship overtakes the waves so 
that the waves appear to actually come from ahead. 

The change in frequency from that of oscillation to wave encounter is analogous 
to the Doppler Shift found in sound waves and electro-magnetic theory. Considering 
ahead seas, the frequency of encounter will be greater than that of the absolute frequency 
and is therefore another critical ingredient to the physics of vessel motion and overall 
understanding of our results. The frequency of encounter and according frequency shifts 
attributed to vessel velocity has a profound influence on ship motions. In particular, 
noting the circular, or concentric, motion amplitudes with high and low, and then high 
again motion amplitude results. These results are due in part to the speed-polar plot data 
representation. With zero velocity at the center, and increasing outwardly in integer 
intervals to a maximum velocity of twenty knots, the frequency shifts are rather evident. 
All statistical assumptions remain the same, as only the point of reference has changed; 
however, the directional spectrum must also be slightly modified to account for such said 
frequency shifts. 

s e = y 2 ix^f = i;f v®.-p¥P*, (45) 

i j 

Realizing the spectrum of encounter to be the limit as d (3 doi tends to zero: 
5[(a>.) (/ ,P;] = ^/8®.8p (46) 

such that: 
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Sco 


1 - 


2co£/ 0 cos P 
8 


Sco (47) 


solving: 

S(co e , p) = S(co,p )/[l-(2cot/ o /g)cos p] (48) 

The null space on the plots is a combination of when a>l, and/or when the 
denominator for the spectrum of encounter is equal to zero (spreading function). 


Table 2. Horizontal Force Comparison 


Horizontal Force Differential (Connection-Ratio) 



5 

0.1 

-28.60% 

0.5 

28.60% 

1.0 

H—' 

5 

0.1 

0.00% 

0.5 

20.00% 

1.0 

sz 

CD 

5 

0.1 

-12.50% 

0.5 

28.60% 

1.0 

'(D 

X 

10 

0.1 

-16.60% 

0.5 

16.60% 

1.0 

CD 

*> 

10 

0.1 

-9.10% 

0.5 

27.30% 

1.0 

> 

10 

0.1 

-14.30% 

0.5 

33.30% 

1.0 

> 

20 

0.1 

-20.00% 

0.5 

20.00% 

1.0 

ci 

o 

20 

0.1 

-12.00% 

0.5 

12.00% 

1.0 

H— 

'c 

20 

0.1 

-14.30% 

0.5 

33.30% 

1.0 

CD 

in 

30 

0.1 

-12.50% 

0.5 

12.50% 

1.0 

30 

0.1 

-18.20% 

0.5 

18.20% 

1.0 


30 

0.1 

-14.30% 

0.5 

14.30% 

1.0 
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V. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

The response of a ship advancing in a stochastic seaway is obviously a 
complicated phenomenon involving the distinct and recognizable interaction between 
hydrodynamic forces and vessel dynamics. Quantifying and measuring either sides of the 
equilibrium equation (39) may never be a perfect science, however, experiment and 
prediction prove theoretical insight is a must for advances in naval architecture. While 
the present discussion involves only distinct evidence for vessel response to waves of 
various, albeit steady, frequencies, we can subsequently still predict the statistics of the 
responses to actual random seaways. 

The measured distribution of hydrodynamic forces along the length of the vessels 
(strip theory predictions) for lower frequencies (longer wavelengths) become greater and 
more adversely affected by the effects of three-dimensionality and forward speed. In all 
cases, lower wave periods have the greatest impact on the vertical and horizontal 
connection forces (on the order of 3.0). For higher frequencies the agreement between 
the strip theory approximation and experimentation are reasonably good for zero forward 
speed, but the agreement between theory and experiment diverges as the speed increases 
(Beck). 

As can be seen by the speed polar plots, forward speed dependence is clearly 
visible in the numerical results. Pitch motion becomes very small in beam seas, but 
become most severe in head seas. The effects of wave frequency and forward speed in 
head seas have been illustrated. In general, a reduction in forward speed will reduce the 
heave and pitch motion in long wavelengths and increase the motion in short waves. The 
motions in short waves (L/L<1.0) are significantly less than in longer wavelengths. 
Vertical connection forces remained relatively the same for increases in the connection 
ratio varying only a maximum of 9.1%. Again, keeping all parameters constant and 
varying the connection ratio, the horizontal components on the other hand ranged from 0- 
33%. The most important results/parameters for seakeeping responses, and amplitude 
motion in random seas are wave period, wave amplitude (significant wave height), and 
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forward speed, and a weaker function of connection ratio for the ranges considered. 
These results compare with results found earlier (Okan), whereby the connection forces 
are not a monotonic function of the towline length (therefore, it is possible to select a 
towline length so that the connection force in minimized). 

For the translational modes, the surge and heave motions were maximized in seas 
forward of the bow. Horizontal plane motions tended to be damped out while vertical 
plane motions grew larger with respect to speed for wave angles forward of beam. For 
sway, it was beam seas that excited the largest motions while yaw was greatest in 
quartering seas. Sway and yaw have no hydrostatic restoring forces displaying high 
amplitude motions in quartering and beam seas. The wave encounter frequencies are also 
much lower in quartering and following seas than in head seas or on the bow and 
evidently lead to larger amplitude sway and yaw motions. 

It must be remembered that the incident and diffracted wave velocity potential’s 
(39) (respective wave velocity potentials before and after contact is made with vessel) 
unsteady components are all higher order and are neglected in the derivation of excitation 
forces using the linear theory assumption. Viscous effects also play a more important 
role in the transverse modes than in the horizontal. For example, if the ship is at an 
apparent angle of attack and coincidently yaws relative to the ship’s forward motion, 
viscous effects will undoubtedly come into effect at the bow. Significant transverse 
rotation, relative to angle of attach, will transpire. However such motions have not been 
incorporated into the ship motion algorithm. Additionally, surge motion is neglected due 
to the simplified head seas case where the vessel surges forward on the surface of the 
wave, also known as surfing. Completing the free-body-diagram on the vessel/wave 
interface, it is easy for one to picture the buoyancy force normal the wave face and 
opposing the vessel body weight component. The included angle between the normal 
force and the weight component is very small, none-the-less the included positive vector 
is a surge force and is clearly present. Sufficiently large tension values, however, may 
stabilize and counter such forces in the coupled two-vessel system 
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Continuing with our discussion, for very long waves the heave motion amplitude 
approaches the wave amplitude and the pitch amplitude is the same as the maximum 
wave slope for such conditions where the vessel moves with the waves (planar motion) 
resulting in lower transverse forces. As the ships are modeled in-tow, they move together 
similarly at the corresponding connection points as a result of the coupling connection. 
As the connection ratio decreases, the leading and trailing ship motions become closely 
matched as the towline length is reduced, while some variation is observable for larger 
towline lengths. Furthermore, the horizontal connection force becomes larger for towing 
lengths (1/L>0.1), while the vertical force decreases (of varying inconsistent magnitude). 
Additionally, there is no clear dependence of the connection force and consequent ship 
motions due to forward speed, there is a large variation with no evident trend. 

The complexities of wave resistance, incident waves and associated singularities 
(wave breaking), as well as some of the nonlinearities associated with the free surface 
problem tend to make the speed polar plots appear more as a glimpse into the amplitude 
motions of six degrees rather than absolute. While the most advanced present-day 
techniques are beyond the capacity of modern-day computers, and will likely be the case 
for years to come, perhaps the simplex method describe here will offer a numerical 
breach into the many unknowns and possibilities of ship response to random impulses. In 
fact, in the open ocean, the evolution and development of waves from ripple to sea state 
ten is likely never to be observed given the obvious explanation that the wave source is, 
in fact, rarely an impulse. Although in a physical sense the ocean wave and its multiple 
sea states cannot be reproduced in a laboratory, from a mathematical standpoint this is 
quite feasible. “In fact, the theoreticians have become so bold as a result of the success 
of their complicated equations that there is danger the study of waves will fall entirely 
into the hands of those who have never [even] been to sea.” (Bascom) 
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B. RECOMMENDATIONS 


The overall effect of vessel control must be realized in order to better determine 
the surge and roll motion amplitudes. Surge is severely affected by propeller thrust, 
rudder design, and steering gear response. Propeller response and the unsteady wake 
behind SLICE would either serve to dampen or excite mode amplitudes for the 
KAIMALINO. Furthermore, the idealized mechanical linkage between the two vessels 
is a major contributor to the coupling mechanics between the vessels and correlates 
directly to roll motion. The damping arising from wave making due to roll is very small 
for this mode; while for heave, yaw, sway, etc. it is a major contributor. This idealized 
linkage assumption, however, should be altered to incorporate a more realistic 
damping/excitation force when considering the roll motion. 

Numerically, the boundary conditions on the hull and on the free surface both 
contained terms based on the steady perturbation potential. In trying to solve this 
boundary value problem, it is determined that there exists an interaction between the 
steady and unsteady components of this potential. Addressing the nonlinearities 
associated with the potential function, in particular, the velocity squared term in the 
pressure/Bernoulli equation and viscous (damping) effects. Similarly, the smooth free 
surface assumption should be reconsidered (incident wave brake). Unique solutions of 
the Laplace equation (Appendix B) require boundary conditions on all surfaces 
surrounding the fluid domain. The major nonlinearities in the general potential flow 
problem are in the free-surface boundary conditions that involve the square of the fluid 
velocities and products of these velocity terms with the unknown free-surface amplitude. 
This potential flow problem still is quite difficult and further underlying assumptions will 
be necessary. For our simplifying assumption of linear theory, the interaction was 
considered negligible and only the free-stream value was deemed necessary. Although 
many in the hydrodynamics field practice this assumption; is this a reasonable hypothesis 
considering our significant wave heights and varying ship velocities? 

Any disturbance of the water surface, including the passage of a ship, produces 
waves. Much of the power expended in propelling the ship goes into such said wave- 
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making processes and while this is a prime concern of the marine engineer in the interest 
of inefficiencies and power losses, it may also be of interest to future proceedings and 
research into using SWATH vessels in-tow. As the ship passes through the water it is 
accompanied by at least three different pressure disturbances on each side, which produce 
several distinct wave trains in consistent order, direction, and symmetry. Lord Kelvin 
investigated the pattern of waves generated by an idealized pressure disturbance 
concentrated at a point and moving in a straight line. The Kelvin wave pattern, as it is 
called, is characterized by (1) diverging waves (a series of curved crests, concave 
outward and lying in echelon position); (2) transverse waves (convex forward and 
perpendicular to the direction of motion); (3) crest intersections (where the diverging and 
transverse waves coincide). It is this trailing pattern that should be further researched and 
incorporated into the two-vessel design. The Volex wave system (ship applications as 
apposed to point disturbance) and pressure deductions should be determined for the 
SLICE twin hull configuration. Do the resultant stern waves contribute to the overall 
wave resistance? When considering the positive pressure wave at the stem, should this 
attribute to added resistance and be incorporated into the connection force? 

Finally, should a more realistic six degree-of-freedom coupling be developed, the 
potential to model this result at the Fleet Numerical Oceanographic Center (FNOC), 
Monterey, CA, in real-time has been addressed to FNOC, and may better ascertain the 
feasibility and validity of the coupling algorithm. (Virtually place the 

SFICE/KAIMAFINO coupling anywhere in the world.) 
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Anecdotal: 


With the advantages of the computational domain come the disadvantages. It has 
been found, without careful observation numerical waves can reflect from the boundary 
conditions and return to and interfere with initial conditions fouling results (Beck and 
Reed, 2001). 

“To impact ship design, computational methods must be integrated into the design 
process. This in turn means that the computer codes must be fast enough that they can 
provide results in no more than an hour or two of computing time. While there has been 
work done on improving the computational efficiency of computer codes for ship 
hydrodynamics problems, a tremendous amount of work still needs to be done. The work 
to date on computational efficiency has resulted in a large enough literature base that this 
subject is worthy of a paper of its own; time and space [and money] obviously preclude 
this.” (Beck and Reed, 2001). 
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APPENDIX A. SPEED POLAR PLOT RESULTS 


Considering the limited memory presently allocated to MS WORD, only a small 
sampling (20 modal solutions) of the 360 member set will be included to illustrate the 
numerical results. If further detail is required of the results outlined in Table 1, please 
email Professor Fotis Papoulias (papoulias@nps.navy.mil) . Mechanical Engineering 
Department, Naval Postgraduate School; Monterey, CA. 
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Vertical Connection Force 


90 


20 



n 0.045 
0.04 

• 10.035 

- • 0.03 

- • 0.025 

- • 0.02 

- -0.015 

■ 0.01 
0.005 


Figure A-l. Vertical Connection Force (HS: 5 feet, T m : 5 sec., 1/L: 0.1) 



Figure A-2. Transverse Connection Force (HS: 5 feet, T m : 5 sec., 1/L: 0.1) 
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Leading Ship, Heave 



Figure A-3. SLICE Heave Amplitude Motion (HS: 5 feet, T m : 5 sec., 1/L: 0.1) 

Trailing Ship, Heave 
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Figure A-4. KAIMALINO Heave Amplitude Motion (HS: 5 feet, T m : 5 sec., 1/L: 

0 . 1 ) 
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Leading Ship, Pitch 



Figure A-5. SLICE Pitch Amplitude Motion (HS: 5 feet, T m : 5 sec., 1/L: 0 

Trailing Ship, Pitch 



Figure A-6. KAIMA L INO Pitch Amplitude Motion (HS: 5 feet, T m : 5 sec., 
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Leading Ship, Sway 



Figure A-7. SLICE Sway Amplitude Motion (HS: 5 feet, T m : 5 sec., 1/L: 0.1) 

Trailing Ship, Sway 
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Figure A-8. KAIMALINO Sway Amplitude Motion (HS: 5 feet, T m : 5 sec., 1/L: 

0 . 1 ) 
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Leading Ship, Yaw 



Figure A-9. SLICE Yaw Amplitude Motion (HS: 5 feet, T m : 5 sec., 1/L: 0.1) 
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Figure A-10. KAIMALINO Yaw Amplitude Motion (HS: 5 feet, T m : 5 sec., 1/L: 

0 . 1 ) 
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Vertical Connection Force 
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Figure A-351. Vertical Connection Force (HS: 30 feet, T m : 20 sec., 1/L: 1.0) 
Transverse Connection Force 

0.14 

0.12 

0.1 

0.08 

0.06 

0.04 

0.02 


Figure A-352. Transverse Connection Force (HS: 30 feet, T m : 20 sec., 1/L: 1.0) 
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Leading Ship, Heave 



Figure A-353. SLICE Heave Amplitude Motion (HS: 30 feet, T m : 20 sec., 1/L: 1.0) 


Trailing Ship, Heave 
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Figure A-354. KAIMALINO Heave Amplitude Motion (HS: 30 feet, T m : 20 sec., 1/L: 

1 . 0 ) 
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Figure A-355. SLICE Pitch Amplitude Motion (HS: 30 feet, T m : 20 sec., 1/L: 1.0) 


Trailing Ship, Pitch 
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Figure A-356. KAIMALINO Pitch Amplitude Motion (HS: 30 feet, T m : 20 sec., 1/L: 

1 . 0 ) 
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Leading Ship, Sway 




Figure A-357. SLICE Sway Amplitude Motion (HS: 30 feet, T m : 20 sec., 1/L: 1.0) 

Trailing Ship, Sway 
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Figure A-358. KA1MA L INO Sway Amplitude Motion (HS: 30 feet, T m : 20 sec., 1/L: 

1 . 0 ) 
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Leading Ship, Yaw 



Figure A-359. SLICE Yaw Amplitude Motion (HS: 30 feet, T m : 20 sec., 1/L: 1.0) 

Trailing Ship, Yaw 



Figure A-360. KAIMALINO Yaw Amplitude Motion (HS: 30 feet, T m : 20 sec., 1/L: 

1 . 0 ) 
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APPENDIX B. FUNDAMENTAL HYDRODYNAMIC THEORY 


The six-degree of freedom motion of the SLICE/KAIMALINO vessel 
configuration in-tow requires a fairly good understanding of the three dimensional 
motion and the corresponding wave hydrodynamics. In order to investigate the motions 
on the ocean surface (particularly wave velocity), we must first develop the basic 
equations (and understanding) of the hydrodynamics of an ideal (nonviscous) fluid. In 
our study, we assumed that the water is incompressible, the density is constant and the 
flow is irrotational. These are standard assumptions in seakeeping studies. We also 
assumed that the ships operate in infinitely deep water. The characterizations are as 
follows: (1) kinematic equations; (2) conservation of mass (continuity equation); (3) 
dynamic equations or Bernoulli’s equation describing irrotational motion of an 
incompressible fluid. 

Using lagrangian techniques, we can describe any number of particles and their 
relative positions with respect to the fluid motion. Letting E = e(v, y,z,t) stand for a 
(scalar) property of the ocean wave medium, the lagrangian rate of change for a particle 
function, E, for some fixed point r G : 

E « = E [xo,y 0 ,Zo,tJ (1) 


When a particle at point r at the time, t, arrives at a neighboring point at the time, 
t+dt, the function E is the following: 


e[x 0 + dx, y„ + d\\f , + dC, ,t a + dt ) = E o + 


"aiU 


dx + 


f dE^ 

\ d yj 0 


d\\f + 


^3E a 
v cte / 


4 + 


^3E a 
\<b a 


dt (2) 


Whereby the time rate of change of E from the lagrangian point of view: 


dE __ dE 
— =v ■ VE+— 
dt dt 


( 3 ) 


where, 
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( 4 ) 


. 3X . 3i|/ , dt, 
v =i^r+j^r+ k ^r 
at at at 

^ .3E . 0E , 3E 

VE = l —+ J — +k — 
ox ay oz 


Thus is described the general kinematic equation of a water particle. 

For our assumptions of an incompressible fluid, which is an accurate model of the 
ocean medium at the surface where the pressure changes are minimal, the continuity 
equation resolves to: 

V • v = 0 (5) 

where v is the displacement velocity which is derivable from the velocity 
potential, (f> , whereby the velocity potential in an incompressible fluid must satisfy 
Laplace’s equation: 

V ■ V4> = 0 (6) 

The surface wave is the manifestation of pressure changes and water-particle 
motions affecting the entire body of fluid- in a column sense all the way through its 
respective depth, surface to ocean floor. The motion of particles under the idealized 
given conditions can be characterized by this velocity potential, <f> . From this function 
all desired wave properties can be determined. 

We shall now derive the Bernoulli Equation by applying Newton’s second law to 
the mass pAx in an element of volume Ax of a fluid. The net force acting on the element 
to give it an acceleration dv/dt includes an external body force, F e (due to gravity), and an 
internal force F (pressure gradients within the fluid). The pressure force between each 
element can be described as follows: 


AF. = -VP Ax (7) 

A F e = - pVflAx (8) 

where Q. is the gravitational potential energy per unit mass, g=-VQ (i.e. positive 
y-axis is directed upward, then fl=gy, and g=-jg. 
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F=ma (9) 


- VP At - pVE>Ax = p —At 

dt 


( 10 ) 


referring to our earlier discussion of equation (3), 


dE „ dE 
— = v ■ VE + — 
dt dt 


( 11 ) 


we may now express the acceleration as follows: 

dv 

- VP- pV£2 = v ■ Vv +— (12) 


When generally discussing waves on the surface of a liquid, it is customary to 
assume that the liquid particle motion is irrotational. This means that the curl of v 
everywhere vanishes, and essentially the mathematical characteristics of eddies, or 
vorticites is lost and obviously not of our concern. When described in an ideal fluid, 
these vortices persist forever to conserve angular momentum; however, in reality vortices 
die out due to fluid viscosity (i.e. smoke rings). Therefore, 

Vxv = 0 (13) 

and, 

(Vxv) xv = v • Vv - Vv ■ v = 0 (14) 

so that: 

V Vv = (Vv)-v = V (y V 2 ) (15) 

The no curl condition is a necessary and sufficient condition that v can be derived from 
the velocity potential (f>. 


v = -V4> (16) 

and since, 


dv 
d t 


dt 


(17) 


equation 10 becomes: 
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VP-VC2 = -V-^-+v( TV 2 ) (18) 
p at 

integration can then be carried through by dotting dr into it and noting that: 
dE = dr ■ VE (19) 

where d E is the difference between two neighboring points separated by dr. 
Equation 18 becomes: 

dP-da = -d^~ + d(jV 2 ) (19) 
p at 

and when integrated gives the relation: 

f dP deb , , 

J—= ^-n-|v 2 (20) 

p at 

and when considered to be incompressible, we resolve to Bernoulli’s Equation: 

p , _ 

-+Q.+ iv 2 = ^- + F(t) (21) 
p at 

whereby the constant of integration, and any time-dependent component has been 

dp 

absorbed into the undetermined —— component. Taking the time derivative and treating 

at 

P as a constant and considering: 

dX dp 


V. 


d t ck 

dv|/ dp 


v = — = - — (22) 
% d t dv ’ 


v_ = 


d c dp 


d t dz 

evaluating y=h, we find: 

dp 1 d 2 p 


V v =-^. - .. ,.2 


(23) 


d y g dV 

showing that velocity depends on position and time, Bernoulli’s equation can now 


be used to compute pressure at any point in the fluid. Solving for the velocity potential 
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function and knowing that it is a function of position and time, we must satisfy the 
Laplace equation: 


3 2 cf) 3 2 cf) 

dx 2 + dy 2 


= 0 (24) 


where; 

^>(x,y,t)=Z(x)-F(0-r(0 (25) 


The spatial functions are determined from separation of variables, and the general 
solution: 


X(x) = Ae iKX + Be ,KX 
Y(y) = Ce Ky +De~ Ky 


(26) 


considering the boundary conditions at y=0: 
deb dY 

v . = - f =- xr ^ =0 <27) 

which requires that C=D, where; 

F(y)= 2Ccosh(Ky) (28) 


similarly the velocity potential has the form: 

() = Acosh k ye iKX T(t) (29) 

and substituting into equation 24, we obtain the equation for T(t): 

d 

——+ (gKtanh k/i)T = 0 (30) 
dt 

which leads to the simple solution of harmonic time dependence with the angular 
frequency: 

co = (gKtanhK/t) 1/2 (31) 
wave velocity is then: 
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—tanh K h 


1/2 


(32) 


for the case of deep water, say when h>?i/2, tanhKh>tanh7t =0.996, so that the 
expression for the wave velocity becomes rather closely: 




, a<h) (33) 
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APPENDIX C. VERTICAL PLANE COUPLING DERIVATION 


Defining the equations of motions for a point on the each of the SLICE and the 
KAIMALINO: 

[A s +A s ]M} + [5 s ]{Ti} + [C s ]{Ti} = {F„} + {/ s } (1) 

[A K + A K M} + [B K M} + [C K ]{r\} = {F ex } + {f K } (2) 

where f s and f K are the connection forces acting on SLICE and KAIMALINO 
respectively, where f K = —f s . The connection forces and absolute motions are 
illustrated in Figure C 1. 



Figure C 1. Connection Force and Motions of the Ships (From: Orhan) 


Continuing with the derivation of Equation (39): 

Ro e 2 (A +A) + m e B + Cm = {F ex } + {/} (3) 

Since the exponential exists in all terms, it is canceled and the equations of 
motion in the frequency domain become as follows: 

A = [-co/(A + A) + zco f 5 + C] (4) 
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We know that A is a 6 x 6 matrix for six degrees of freedom system: 

Ap = {F} (5) 

rf = inv{A){F} (6) 

where; 

{F} = {F a } + {f} (7) 

Ship motions due to regular waves of a given wavelength and direction are now 
determined for a given forward speed (V). Motions in vertical and horizontal planes are 
decoupled for our discussion and may be solved as two distinct 3x3 systems. The heave 
and pitch modes and the resultant interactions will be derived from this point forward 
(sway and yaw are similar). The expanded equations of motion in two degrees of 
freedom become as: 

A 33 Ajj 

_Ai3 ^55 

Now we can define the equations of motion in vertical plane (two degrees of 
freedom) for the two vessels as: 

A33..T1 3 S . + A35.sr| 5iS , = F 3S + f s (12) 

A53,jT ) 3 s + A 55.51) 5 5 — F 5 S ~ f s X s (13) 

A33.kT\ 3 k + A35,kT\ 5 k = F 3 k + f K (14) 

A53.a:T1 3 k + A55,/rT ] 5 K — F 5 K — f s X K (15) 


fn,- 


' F 3 +f 



F 5 + fx s _ 


( 11 ) 
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where x is the distance between the connection point and the center of gravity of 
that vehicle. In order to solve these four equations, we make the following substitution: 


0 3,S = 

M 3 .S 

+v 

3 ,sfs 

His 

Ms.s 

+v, 

5 ,sfs 

n - _ 

M 3 , A 

+v 

3, A f K 

TIs.a = 

M 5 , A 

+v 

5, A f A 

where; 





(16) 


fs=-f K = -f (17) 

Using Cramer’s Rule, the absolute motions can be described by (Nash): 
^ = r \ 3 ,s 

- - U°) 

J A _r l3,A' ~^i,K X K 

When we combine (13) and (14) as: 

^>S = M3,S 3,S f ~ ^5,S X S +V 5 ,S X sf 
^3 A = M 3 ,A 3 ,k/ — \^5,K X K ~ V 5 ,K X K f 

Let a in terms of p ’s and b in terms of v’s be defined as: 

a — M3,5 _ M3.A Ms.A’^'A ~ M'S, S X S 
3,S ~^ 3,K ^5,K X K 5,S X S 


( 20 ) 


Therefore, our equation of motion becomes: 
£>s ~^>k = a~bf (21) 

Again referring to Figure 2, 

/ = T ^s~^K ( 22 ) 
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Combining (17) and (18) gives us our final equation, which gives us the 
connection force between two ships: 


f=-r— ( 23 ) 

—+b 

T 



Figure C 2. KAIMALINO Tension vs. Speed Plot (After: Lockheed Martin) 
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APPENDIX D. MATLAB COUPLING ALGORITHM 


% Horizontal Plane 
% Dimensional version (U.S. units) 

% Contour plots (heading/speed) 

% Two parameter Bretschneider spectrum - Short crested seas 


% 

% Get run info 

% 

HS=input('Significant Wave Height (feet) = '); 

T_m=input('Modal Period 
l=input('Length (l/L) = '); 

% 

omega_m=2*pi/T_m; 

HS_string=num2str(HS); 
warning off 

% 

lambda_min=20; 
lambda_max=1000; 
delta_lambda=20; 
rho=1.9905; 
zeta=l; 

L=105; 
g=32.2; 
x_s=-46; 
x_k=+40; 
beta_incr=5; 
motion_ratio=2; 

% 

% The matdata output files default to the vertical only format when the 
% heading angle is 0 or 180 degrees. 

% Set up file reading format. 

% 

% 

% GENERAL DATA 

% 

lambda=lambda_min:delta_lambda:lambda_max; % Vector of wavelengths 

wavenumber=2.0*pi./lambda; % Wave number 

omega=sqrt(wavenumber*g); %Wave frequency 

period=2.0*pi./omega; 

omega=omega'; 

filesize=size( lambda); 

if beta_incr= = 15 

load data_ 15; % Load main data file, 15 deg. heading increments 

beta_incr_values=13; 
end 

if beta_incr==5 

load data_5; % Load main data file, 5 deg. heading increments 

beta_incr_values=37; 
end 

% 


% Min wave length (ft) 

% Max wave length (ft) 

% Wave length increment (ft) 

% Water density 
% Regular wave height 
% Reference length for nondimensionalization 
% Gravitational constant 
% FRONT SHIP attachment point 
% REAR SHIP attachment point 
% Increment in sea direction (deg) 

% Ratio for motion plotting comparison 
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% Random wave calculations 

% Bretschneider spectrum - weighted by the spreading function 

% 

iSpeed=0; 

for V_knots=0:l:20, % Loop on speed 

iSpeed=iSpeed + l; 

% 

% Get tension from curvefitting data. 

% Applicable for speeds between 1 and 20 ft/sec. 

% 

T =-1.762* V_knots ~4+63.675* V_ knots ~ 3-580.8*V_ knots''2+ 2485.9* V_knots-34.047; 

% 

V_string=num2str(V_knots) 

V=V_knots* 1.6878; % Convert to ft/sec 

A=(1.25/4)*(omega_m /v 4)*(HS / '2); 

B=1.25*omega_m'M; 

S_main=(A./omega.''5).*exp(- B./omega.'M); 

% 

ibeta=0; 

for beta=0:beta_incr:360, % Loop on sea direction 

beta_set=beta; 
if beta>180 

beta_set=360-beta; 

end 

ibeta=ibeta + l; 
beta; 

% 

% Spreading between (beta-90) and (beta+90) 

% 

beta_spread_low =beta_set - 90; 
beta_spread_high=beta_set + 90; 
ibeta_spread=0; 

% 

% Loop between (beta-90) and (beta + 90) 

% Use either 15 deg increments (13 values total) or 
% 5 deg increments (total of 37 values) 

% 

for beta_spread=beta_spread_low:beta_incr:beta_spread_high, 
ibeta_spread=ibeta_spread + 1; 
beta_spread_read = beta_spread; 
beta_spread_vector(ibeta_spread)=beta_spread; 
if beta_spread<0 

beta_spread_read=-beta_spread; 

end 

if beta_spread>180 

beta_spread_read=360-beta_spread; 

end 

beta_string=num2str(beta_spread_read); 

trigg=30; 

f2loc=26; f6loc=30; 
if beta_spread_read ==0 
trigg=27; 

f2loc=25; f6loc=27; 
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elseif beta_spread_read = = 180 
trigg=27; 

f2loc=25; f6loc=27; 
end 

lambda_size=trigg*filesize(2); 

% 

% Load FRONT SHIP data file msvhV_beta.txt 

% 

load_filename=strcat('msvh',V_string,'_',beta_string); 

filename_s=eval(load_filename); 

% 

% Load REAR SHIP data file 

% 

load_filename=strcat( , mkvh',V_string,'_',beta_string); 

filename_k=eval(load_filename); 

omegae=omega'-wavenumber*V*cos(beta_spread*pi/180); %Freq. of Encounter 

omegae=omegae'; 

periode=2.0*pi./omegae; 

% 

% HORIZONTAL PLANE RESPONSE CALCULATIONS 

% 

% SLICE 

% 

% Set mass matrix elements 

% 

M22s=filename_s(2:trigg:lambda_size,2); 

M26s=filename_s(2:trigg:lambda_size,6); 

M62s=filename_s(6:trigg:lambda_size,2); 

M66s=filename_s(6:trigg:lambda_size,6); 

% 

% Added mass terms 

% 

A22s=filename_s(8:trigg:lambda_size,2); 

A26s=filename_s(8:trigg:lambda_size,6); 

A62s=filename_s(12:trigg:lambda_size,2); 

A66s=filename_s(12:trigg:lambda_size,6); 

% 

% Damping terms 

% 

B22s=filename_s(14:trigg:lambda_size,2); 

B26s=filename_s(14:trigg:lambda_size,6); 

B62s=filename_s(18:trigg:lambda_size,2); 

B66s=fiiename_s(18:trigg:lambda_size,6); 

% 

% Hydrostatic terms 

% 

C22s=filename_s(20:trigg:lambda_size,2); 

C26s=fiiename_s(20:trigg:lambda_size,6); 

C62s=filename_s(24:trigg:lambda_size,2); 

C66s=filename_s(24:trigg:lambda_size,6); 

% 

if beta_spread ==0 

F2s_t=zeros(50,l); F6s_t=zeros(50,l); 
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elseif beta_spread = = 180 

F2s_t=zeros(50,l); F6s_t=zeros(50,l); 
else 

% 

% Total exciting forces 

% 

F2s_t_amp=filename_s(f2loc:trigg:lambda_size,5); 

F6s_t_amp=filename_s(f6loc:trigg:lambda_size,5); 

F2s_t_pha=filename_s(f2loc:trigg:lambda_size,6); 

F6s_t_pha=filename_s(f6loc:trigg:lambda_size,6); 

F2s_t=F2s_t_amp.*exp(i*F2s_t_pha* pi/180.0); 

F6s_t=F6s_t_amp.*exp(i*F6s_t_pha.*pi/180.0); 

% 

% Froude/Krylov exciting forces 

% 

F2s_f_amp=filename_s(f2loc:trigg:lambda_size,l); 

F6s_f_amp=filename_s(f6loc:trigg:lambda_size,l); 

F2s_f_pha=filename_s(f2loc:trigg:lambda_size,2); 

F6s_f_pha=fiiename_s(f6loc:trigg:lambda_size,2); 

F2s_f=F2s_f_amp.*exp(i*F2s_f_pha.*pi/180.0); 

F6s_f=F6s_f_amp.*exp(i*F6s_f_pha* pi/180.0); 

% 

% Diffraction exciting forces 

% 

F2s_d_amp=filename_s(f2loc:trigg:lambda_size,3); 

F6s_d_amp=filename_s(f6loc:trigg:lambda_size,3); 

F2s_d_pha=filename_s(f2loc:trigg:lambda_size,4); 

F6s_d_pha=filename_s(f6loc:trigg:lambda_size,4); 

F2s_d=F2s_d_amp.*exp(i*F2s_d_pha.* pi/180.0); 

F6s_d=F6s_d_amp.*exp(i*F6s_d_pha.*pi/180.0); 

% 

end 

% KAI MALI NO 

% 

% Set mass matrix elements 

% 

M22k=filename_k(2:trigg:lambda_size,2); 

M26k=filename_k(2:trigg:lambda_size,6); 

M62k=filename_k(6:trigg:lambda_size,2); 

M66k=filename_k(6:trigg:lambda_size,6); 

% 

% Added mass terms 

% 

A22k=filename_k(8:trigg:lambda_size,2); 

A26k=filename_k(8:trigg:lambda_size,6); 

A62k=filename_k(12:trigg:lambda_size,2); 

A66k=filename_k(12:trigg:lambda_size,6); 

% 

% Damping terms 

% 

B22k=filename_k(14:trigg:lambda_size,2); 

B26k=filename_k(14:trigg:lambda_size,6); 

B62k=filename_k(18:trigg:iambda_size,2); 
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B66k=filename_k(18:trigg:lambda_size,6); 

% 

% Hydrostatic terms 

% 

C22k=filename_k(20:trigg:lambda_size,2); 
C26k=filename_k(20:trigg:lambda_size,6); 
C62k=filename_k(24:trigg:lambda_size,2); 
C66k=filename_k(24:trigg:lambda_size,6); 
if beta_spread ==0 

F2k_t=zeros(50,l); F6k_t=zeros(50,l); 
elseif beta_spread = = 180 

F2k_t=zeros(50,l); F6k_t=zeros(50,l); 
else 

% 

% Total exciting forces 

% 

F2k_t_amp=filename_k(f2loc:trigg:lambda_size,5); 
F6k_t_amp=filename_k(f6loc:trigg:lambda_size,5); 
F2k_t_pha=filename_k(f2loc:trigg:lambda_size,6); 
F6k_t_pha=filename_k(f6ioc:trigg:lambda_size,6); 
F2k_t=F2k_t_amp.*exp(i*F2k_t_pha.* pi/180.0); 
F6k_t=F6k_t_amp.*exp(i*F6k_t_pha.* pi/180.0); 

% 

% Froude/Krylov exciting forces 

% 

F2k_f_amp=filename_k(f2loc:trigg:lambda_size,l); 
F6k_f_amp=filename_k(f6loc:trigg:lambda_size,l); 
F2k_f_pha=filename_k(f2loc:trigg:lambda_size,2); 
F6k_f_pha=filename_k(f6loc:trigg:lambda_size,2); 
F2k_f=F2k_f_amp.*exp(i*F2k_f_pha.* pi/180.0); 
F6k_f=F6k_f_amp.*exp(i*F6k_f_pha.* pi/180.0); 

% 

% Diffraction exciting forces 

% 

F2k_d_amp=filename_k(f2loc:trigg:lambda_size,3); 
F6k_d_amp=filename_k(f6loc:trigg:lambda_size,3); 
F2k_d_pha=filename_k(f2loc:trigg:lambda_size,4); 
F6k_d_pha=filename_k(f6loc:trigg:lambda_size,4); 
F2k_d=F2k_d_amp.*exp(i*F2k_d_pha.* pi/180.0); 
F6k_d=F6k_d_amp.*exp(i*F6k_d_pha.* pi/180.0); 
end 
% 

% MATCHING CONDITION 

% 

A22bar_s=-(omegae.~2).*(M22s+A22s)+i*omegae.*B22s+C22s; 
A26bar_s=-(omegae.~2).*(M26s+A26s)+i*omegae.*B26s+C26s; 
A62bar_s=-(omegae.~2).*(M62s+A62s)+i*omegae.*B62s+C62s; 
A66bar_s=-(omegae.~2).*(M66s+A66s)+i*omegae.*B66s+C66s; 
A22bar_k=-(omegae.''2).*(M22k+A22k) + i*omegae.*B22k+C22k 
A26bar_k=-(omegae.~2).*(M26k+A26k)+i*omegae.*B26k+C26k 
A62bar_k=-(omegae.~2).*(M62k+A62k)+i*omegae.*B62k+C62k 
A66bar_k=-(omegae.~2).*(M66k+A66k)+i*omegae.*B66k+C66k 
% 
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A26_s_den=(A22bar_s.*A66bar_s-A26bar_s.*A62bar_s); 

A26_k_den = (A22bar_k.*A66bar_k-A62bar_k.*A26bar_k); 

% 

mu2_s=(A66bar_s.*F2s_t-A26bar_s.*F6s_t)./A26_s_den; 

nu2_s=(A66bar_s-A26bar_s*x_s)./A26_s_den; 

mu6_s=(A22bar_s.*F6s_t-A62bar_s.*F2s_t)./A26_s_den; 

nu6_s=(A22bar_s*x_s-A62bar_s)./A26_s_den; 

mu2_k=(A66bar_k.*F2k_t-A26bar_k.*F6k_t)./A26_k_den; 

nu2_k=(A66bar_k-A26bar_k*x_k)./A26_k_den; 

mu6_k=(A22bar_k.*F6k_t-A62bar_k.*F2k_t)./A26_k_den; 

nu6_k=(A22bar_k*x_k-A62bar_k)./A26_k_den; 

% 

a=mu2_s+mu6_s*x_s-mu2_k-mu6_k*x_k; 
b=nu2_s+nu6_s*x_s+nu2_k+nu6_k*x_k; 
f=a./(I/T + b); 

% 

f_s=-f; % Connection force on SLICE 

f_k=f; % Connection force on KAIMALINO 

eta2_s=mu2_s+nu2_s.*f_s; % SLICE sway 

eta6_s=mu6_s+nu6_s.*f_s; % SLICE yaw 

eta2_k=mu2_k+nu2_k.*f_k; % KAI MALI NO sway 

eta6_k=mu6_k+nu6_k.*f_k; % KAIMALINO yaw 

xi_s=eta2_s+eta6_s*x_s; % SLICE motion at connection 

xi_k=eta2_k+eta6_k*x_k; % KAI MALI NO motion at connection 

xi0_s=mu2_s+mu6_s*x_s; % SLICE motion at connection for zero f 

xi0_k=mu6_k+mu6_k*x_k; % KAIMALINO motion at connection for zero f 

delta_beta(ibeta_spread)=beta_set-beta_spread; 

S=S_main *(2/pi)*(cos((beta_set-beta_spread)*pi/180))''2; % Mult.by (2/pi)cos^2 

Se=S./abs((l-(2.0/g)*omega*V*cos(beta_spread*pi/180))); % S(w) to S(we) 

% 

% Define response spectra and store as functions of theta 

% 

Sf(:,ibeta_spread)=((abs(f)).^2).*Se; 

Sxi_s(:,ibeta_spread) = ((abs(xi_s)).^2).*Se; 
Sxi_k(:,ibeta_spread)=((abs(xi_k)).~2).*Se; 
SxiO_s(:,ibeta_spread)=((abs(xiO_s)).^2).*Se; 

SxiO_k(:,ibeta_spread) = ((abs(xiO_k)).''2).*Se; 

SF2s_t(:,ibeta_spread) = ((abs(F2s_t)).^2).*Se; 

SF2k_t(:,ibeta_spread) = ((abs(F2k_t)). /v 2).*Se; 

Seta2_s(:,ibeta_spread) = ((abs(eta2_s)).^2).*Se; 
Smu2_s(:,ibeta_spread)=((abs(mu2_s)).^2).*Se; 

Seta6_s(:,ibeta_spread) = ((abs(eta6_s)).^2).*Se; 
Smu6_s(:,ibeta_spread)=((abs(mu6_s)). ~2).*Se; 
Seta2_k(:,ibeta_spread)=((abs(eta2_k)).^2).*Se; 

Smu2_k(:,ibeta_spread) = ((abs(mu2_k)).^2).*Se; 
Seta6_k(:,ibeta_spread)=((abs(eta6_k)).^2).*Se; 

Smu6_k(:,ibeta_spread) = ((abs(mu6_k)).^2).*Se; 
end % Loop on beta_spreading ends 

% 

% Initialize variables before integrating 

% 

Sf_i=0; 

Sxi_s_i=0; 
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Sxi_k_i=0; 

Sxi0_s_i=0; 

Sxi0_k_i=0; 

SF2s_t_i = 0; 

SF2k_t_i=0; 

Seta 2_s_i=0; 

Smu2_s_i=0; 

Seta6_s_i=0; 

Smu6_s_i=0; 

Seta2_k_i=0; 

Smu2_k_i=0; 

Seta6_k_i=0; 

Smu6_k_i=0; 

% 

% Integrate response spectra over omega and theta 

% 

for J=2:l:beta_incr_values, % 13 for 15 deg increments, 37 for 5 deg increments 
for I =2:l:filesize(2), 

% 

delta_omega=abs(omegae(l-l)-omegae(l)); 

delta_theta=beta_incr*pi/180; 

% 

Sf_sum=Sf(l ,J )+Sf(l-l,J)+Sf(I ,J-1)+Sf(I- 1,J-1); 

Sxi_s_sum=Sxi_s(l,J )+Sxi_s(l-l,J )+Sxi_s(l,J -l)+Sxi_s(l-l,J-1); 

Sxi_k_sum=Sxi_k( I ,J)+Sxi_ k( I-1 ,J) + Sxi_ k( I ,J -1) + Sxi_ k( I-1 ,J -1); 

SxiO s sum=SxiO_s(I ,J) + SxiO_s(I-1 ,J)+ SxiO_s(l,J-1) + SxiO_s( 1-1 ,J-1); 
SxiO_k_sum=SxiO_k(l ,J) + SxiO_k(l-l,J) + SxiO_k(l J-1) + SxiO_k(l-l,J-l); 
SF2s_t_sum=SF2s_t(l ,J) + SF2s_t(l-l,J) + SF2s_t(l ,J-1) + SF2s_t(l-1J-1); 
SF2k_t_sum=SF2k_t(l,J) + SF2k_t(l-l,J) + SF2k_t(l,J-l) + SF2k_t( 1-1 ,J-1); 
Seta2_s_sum=Seta2_s(l,J) + Seta2_s(l-l,J) + Seta2_s(l ,J -1) + Seta2_s( I-1 ,J -1); 
Smu2_s_sum=Smu2_s(l,J) + Smu2_s(l-l,J) + Smu2_s(l,J-l) + Smu2_s(l-l,J-1); 
Seta6_s_sum=Seta6_s(l ,J) + Seta6_s(l-l,J) + Seta6_s(l ,J -1) + Seta6_s(I-1 ,J-1); 
Smu6_s_sum=Smu6_s(l,J) + Smu6_s(l-l,J) + Smu6_s(l,J-l) + Smu6_s(l-l,J-1); 
Seta2_k_sum=Seta2_k(l,J) + Seta2_k(l-l,J) + Seta2_k(l,J-l) + Seta2_k(l-l,J-1); 
Smu2_k_sum=Smu2_k(l,J) + Smu2_k( 1-1 ,J) + Smu2_k(l,J-l) + Smu2_k(l-l,J-l); 
Seta6_k_sum=Seta6_k(l,J) + Seta6_k(l-l,J) + Seta6_k(l,J-l) + Seta6_k(l-l,J-l); 
Smu6_k_sum=Smu6_k(l ,J) + Smu6_k(l-l,J) + Smu6_k(l,J-l) + Sm u6_ k( I-1 ,J -1); 
% 

Sf_i=Sf_i + 0.25*delta_omega*delta_theta*Sf_sum; 

Sxi_s_i=Sxi_s_i + 0.25*delta_omega*delta_theta*Sxi_s_sum; 

Sxi_k_i=Sxi_k_i + 0.25*delta_omega*delta_theta*Sxi_k_sum; 

SxiO_s_i=SxiO_s_i + 0.25*delta_omega*delta_theta*Sxi0_s_sum; 
SxiO_k_i=SxiO_k_i + 0.25*delta_omega*delta_theta*Sxi0_k_sum; 
SF2s_t_i=SF2s_t_i + 0.25*delta_omega*delta_theta*SF2s_t_sum; 
SF2k_t_i=SF2k_t_i + 0.25*delta_omega*delta_theta*SF2k_t_sum; 
Seta2_s_i=Seta2_s_i + 0.25*delta_omega*delta_theta*Seta2_s_sum; 
Smu2_s_i=Smu2_sJ + 0.25*delta_omega*delta_theta*Smu2_s_sum; 

Seta 6_s_i=Seta 6_s_i + 0.25*delta_omega*delta_theta*Seta6_s_sum; 
Smu6_s_i=Smu6_s_i + 0.25*delta_omega*delta_theta*Smu6_s_sum; 
Seta2_k_i=Seta2_k_i + 0.25*deita_omega*deita_theta*Seta2_k_sum; 

Smu2_k_l =Smu2_k_i + 0.25*delta_omega*deita_theta*Smu2_k_sum; 
Seta6_k_i=Seta6_k_i + 0.25*delta_omega*delta_theta*Seta6_k_sum; 
Smu6_k_i=Smu6_k_i + 0.25*delta_omega*delta_theta*Smu6_k_sum; 
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end 

end 

% 

% RMS values 

% 

RMS_f = sqrt(Sf_i); 

RMS_xi_s = sqrt(Sxi_s_i); 

RMS_xi_k = sqrt(Sxi_k_i); 

RMS_xiO_s = sqrt(SxiO_s_i); 

RMS_xiO_k = sqrt(SxiOki); 

RMS_F2s_t = sqrt(SF2s_t_i); 

RMS_F2k_t = sqrt(SF2k_t_i); 

RMS_eta2_s= sqrt(Seta2_s_i); 

RMS_mu2_s = sqrt(Smu2_s_i); 

RMS_eta6_s= sqrt(Seta6_s_i); 

RMS_mu6_s = sqrt(Smu6_s_i); 

RMS_eta2_k= sqrt(Seta2_k_i); 

RMS_mu2_k = sqrt(Smu2_kJ); 

RMS_eta6_k= sq rt( Seta 6_ k_ i); 

RMS_mu6_k = sqrt(Smu6_k_i); 

% 

RMS_f_vector(ibeta,iSpeed) = RMS_f/(rho*g*L''2); 
RMS_eta2_s_vector(ibeta,iSpeed) = RMS_eta2_s/RMS_mu2_s; 
RMS_eta6_s_vector(ibeta,iSpeed)= RMS_eta6_s/RMS_mu6_s; 
RMS_eta2_k_vector(ibeta,iSpeed) = RMS_eta2_k/RMS_mu2_k; 
RMS_eta6_k_vector(ibeta,iSpeed) = RMS_eta6_k/RMS_mu6_k; 

% 

if RMS_eta2_s_vector(ibeta,iSpeed) > motion_ratio 
RMS_eta2_s_vector(ibeta,iSpeed) = motion_ratio; 
end 

if RMS_eta6_s_vector(ibeta,iSpeed) > motion_ratio 
RMS_eta6_s_vector(ibeta,iSpeed) = motion_ ratio; 
end 

if RMS_eta2_k_vector(ibeta,iSpeed) > motion_ratio 
RMS_eta2_k_vector(ibeta,iSpeed) = motion_ratio; 
end 

if RMS_eta6_k_vector(ibeta,iSpeed) > motion_ratio 
RMS_eta6_k_vector(ibeta,iSpeed) = motion_ratio; 
end 

end % Loop on beta ends 
end % Loop on speed ends 

% 

% Save results 

% 

clear msv* 
clear mkv* 

I_string=num2str(l); 

T_m_string=num2str(T_m); 

save_filename=strcat( l h_speed_Bretch_ l ,HS_string, l _',T_m_string, l _ l ,l_string, , .mat'); 

save(save_filename); 
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APPENDIX E. MATLAB AMPLITUDE PLOT ALGORITHM 


% 

% Contour plots - Bretchneider short crested seas 

% 

type=input('Plot type (l=Speed/heading) (2=Waveheight/heading) = '); 

% 

% Speed-Heading plots 

% 

if type= = l 

HS =input('Significant Wave Height (feet) = '); 

T_m =input('Modal Period (sec) = '); 

I =input('Length (l/L) = '); 

HS_string=num2str(HS); 

I_string=num2str(l); 

T_m_string=num2str(T_m); 

% 

% Horizontal plane 

% 

loadf ilename=strcat( 'h_speed_Bretch_HS_string,'_ ',T_ mstring,l_string,'. mat 1 ); 
load(load_filename); 

% 

figure(l) 

[th,r]=meshgrid(( 0: beta_incr: 360)* pi/180,0:1:20); 

[X,Y]=pol2cart(th,r); 
h=polar(th,r);delete(h); 
hold on 

c_p=[0.000:0.005:1.0]; 
contour(X',Y',RMS_f_vector,c_p),colorbar 
title(Transverse Connection Force') 

% 

figure(2) 

[th,r]=meshgrid((0: beta_incr: 360)* pi/180,0:1:20); 

[X,Y]=pol2cart(th,r); 
h=polar(th,r);delete(h); 
hold on 

c_p=[ 0.000:0.010: motion_ratio]; 

contour(X',Y',RMS_eta2_s_vector,c_p),caxis([0 motion_ratio]),colorbar 
title('Leading Ship, Sway') 

% 

figure(3) 

[ th,r]=meshgrid((0: beta_incr: 360)*pi/180,0:l:20); 

[X,Y]=pol2cart(th,r); 
h=polar(th,r);delete(h); 
hold on 

c_p=[ 0.000:0.010: motion_ratio]; 

contour(X',Y',RMS_eta6_s_vector,c_p),caxis([0 motion_ratio]),colorbar 
title('Leading Ship, Yaw') 

% 

figure(4) 

[ th,r]=meshgrid( (0: betaincr: 360)*pi/180,0:l:20); 

[X,Y]=pol2cart(th,r); 
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h=polar(th,r);delete(h); 
hold on 

c_p=[ 0.000:0.010: motion_ ratio ]; 

contour(X',Y',RMS_eta2_k_vector,c_p),caxis([0 motion_ratio]),colorbar 
title('Trailing Ship, Sway 1 ) 

% 

figure(5) 

[th,r]=meshgrid((0:betajncr:360)*pi/180,0:l:20); 

[X,Y]=pol2cart(th,r); 
h=polar(th,r);delete(h); 
hold on 

c_p=[ 0.000:0.010: motion_ ratio ]; 

contour(X',Y',RMS_eta6_k_vector,c_p),caxis([0 motion_ratio]),colorbar 
title('Trailing Ship, Yaw') 

% 

% Vertical plane 

% 

load_filename=strcat('v_speed_Bretch_',HS_string,'_',T_m_string,'_',l_string,'.mat'); 

load(load_filename); 

% 

figure(6) 

[th,r]=meshgrid( (0: betaincr: 360)* pi/180,0:1:20); 

[X,Y]=pol2cart(th,r); 
h=polar(th,r);delete(h); 
hold on 

c_p=[0.000:0.005:1.0]; 
contour(X',Y',RMS_f_vector,c_p),colorbar 
title]'Vertical Connection Force 1 ) 

% 

figure] 7) 

[th,r]=meshgrid((0: beta_incr: 360)* pi/180,0:1:20); 

[X,Y]=pol2cart(th,r); 
h=polar(th,r);delete(h); 
hold on 

c_p=[ 0.000:0.010: motion_ ratio ]; 

contour(X',Y',RMS_eta3_s_vector,c_p),caxis([0 motion_ratio]),colorbar 
title]'Leading Ship, Heave') 

% 

figure] 8) 

[th,r]=meshgrid( (0: beta_incr: 360)* pi/180,0:1:20); 

[X,Y]=pol2cart(th,r); 
h=polar(th,r);delete(h); 
hold on 

c_p=[0.000:0.010: motion_ ratio ]; 

contour(X',Y',RMS_eta5_s_vector,c_p),caxis([0 motion_ratio]),colorbar 
title('Leading Ship, Pitch') 

% 

figure(9) 

[th,r]=meshgrid((0: beta_incr: 360)* pi/180,0:1:20); 

[X,Y]=pol2cart(th,r); 
h=polar(th,r);delete(h); 
hold on 

c_p=[0.000:0.010: motionratio]; 
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contour(X',Y',RMS_eta3_k_vector,c_p),caxis([0 motion_ratio]),colorbar 
title('Trailing Ship, Heave') 

% 

figure(lO) 

[th,r]=meshgrid((0: beta_incr: 360)* pi/180,0:1:20); 

[X,Y]=pol2cart(th,r); 
h=polar(th,r);delete(h); 
hold on 

c_p=[0.000:0.010: motion_ratio]; 

contour(X',Y',RMS_eta5_k_vector,c_p),caxis([0 motion_ratio]),colorbar 
title('Trailing Ship, Pitch 1 ) 
end 

if type==2 

V =input('Speed (knots) = '); 

T_m =input('Modal Period (sec) = '); 

I =input('Length (l/L) = '); 

V_string=num2str(V); 

I_string=num2str(l); 

T_m_string=num2str(T_m); 

% 

% Horizontal plane 

% 

load_filename=strcat('h_waveheight_Bretch_',V_string,'_',T_m_string,'_',l_string,'.mat'); 

load(load_filename); 

% 

figure(l) 

[th,r]=meshgrid((0:beta_incr:360)*pi/180,0.5:0.5:30); 

[X,Y]=pol2cart(th,r); 
h=polar(th,r);delete(h); 
hold on 

c_p=[0.000:0.005:1.0]; 
contour(X',Y',RMS_f_vector,c_p),colorbar 
title(Transverse Connection Force') 

% 

figure(2) 

[th, r]=meshgrid( (0: beta_incr: 360)*pi/180,0.5:0.5:30); 

[X,Y]=pol2cart(th,r); 
h=polar(th,r);delete(h); 
hold on 

c_p=[0.000:0.010: motion_ratio]; 

contour(X',Y',RMS_eta2_s_vector,c_p),caxis([0 motion_ratio]),colorbar 
title('Leading Ship, Sway') 

% 

figure(3) 

[th,r]=meshgrid( (0: beta_incr: 360)* pi/180,0.5:0.5:30); 

[X, Y]=pol2cart(th, r); 
h=polar(th,r);delete(h); 
hold on 

c_p=[ 0.000:0.010: motion_ratio]; 

contour(X',Y',RMS_eta6_s_vector,c_p),caxis([0 motion_ratio]),colorbar 
title('Leading Ship, Yaw') 

% 

figure(4) 
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[th,r]=meshgrid((0: betaincr: 360)*pi/180,0.5:0.5:30); 

[X,Y]=pol2cart(th,r); 
h=polar(th,r);delete(h); 
hold on 

c_p=[ 0.000:0.010: motion_ ratio ]; 

contour(X',Y',RMS_eta2_k_vector,c_p),caxis([0 motion_ratio]),colorbar 
title('Trailing Ship, Sway 1 ) 

% 

figure! 5) 

[th,r]=meshgrid((0:beta_incr:360)*pi/180,0.5:0.5:30); 

[X,Y]=pol2cart(th,r); 
h=polar(th,r);delete(h); 
hold on 

c_p=[0.000:0.010: motion_ ratio ]; 

contour(X',Y',RMS_eta6_k_vector,c_p),caxis([0 motion_ratio]),colorbar 
title!'Trailing Ship, Yaw') 

% 

% Vertical plane 

% 

load_filename=strcat( l v_waveheight_Bretch_ l ,V_string,'_',T_m_string,'_ l ,l_string, , .mat l ); 

load(load_filename); 

% 

figure! 6) 

[th,r]=meshgrid( (0: beta_ incr: 360)*pi/180,0.5:0.5:30); 

[X,Y]=pol2cart(th,r); 
h=polar(th,r);deiete(h); 
hold on 

c_p=[0.000:0.005:1.0]; 
contour(X',Y',RMS_f_vector,c_p),colorbar 
title!'Vertical Connection Force') 

% 

figure! 7) 

[th,r]=meshgrid( (0: beta_ incr: 360)*pi/180,0.5:0.5:30); 

[X,Y]=pol2cart(th,r); 
h=polar(th,r);deiete(h); 
hold on 

c_p=[0.000:0.010: motionratio]; 

contour(X',Y',RMS_eta3_s_vector,c_p),caxis([0 motion_ratio]),colorbar 
title!'Leading Ship, Heave') 

% 

figure! 8) 

[th,r]=meshgrid((0:beta_incr:360)*pi/180,0.5:0.5:30); 

[X,Y]=pol2cart(th,r); 
h=polar(th,r);delete(h); 
hold on 

c_p=[0.000:0.010: motion_ ratio ]; 

contour(X',Y',RMS_eta5_s_vector,c_p),caxis([0 motion_ratio]),colorbar 
title('Leading Ship, Pitch') 

% 

figure(9) 

[th,r]=meshgrid( (0: beta_ incr: 360)*pi/180,0.5:0.5:30); 

[X,Y]=pol2cart(th,r); 

h=polar(th,r);delete(h); 


62 



hold on 

c_p=[0.000:0.010: motion_ratio]; 

contour(X',Y',RMS_eta2_k_vector,c_p),caxis([0 motion_ratio]),colorbar 
title(Trailing Ship, Heave') 

% 

figure(lO) 

[th,r]=meshgrid((0:beta_incr:360)*pi/180,0.5:0.5:30); 

[X,Y]=pol2cart(th,r); 

h=polar(th,r);delete(h); 

hold on 

c_p=[0.000:0.010: motion_ratio]; 

contour(X',Y',RMS_eta5_k_vector,c_p),caxis([0 motion_ratio]),colorbar 
title('Trailing Ship, Pitch 1 ) 
end 
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